Tunneling- induced damping of phase coherence revivals in deep optical lattices 
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We consider phase coherence collapse and revival in deep optical lattices, and calculate within the 
Bose-Hubbard model the revival amplitude damping incurred by a finite tunneling coupling of the 
lattice wells (after sweeping from the superfluid to the Mott phase). Deriving scaling laws for the 
corresponding decay of first-order coherence revival in terms of filling factor, final lattice depth, and 
number of tunneling coupling partners, we estimate whether revival-damping related to tunneling 
between sites can be or even has already been observed in experiment. 
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Introduction. The idea of implementing the Bose- 
Hubbard model with cold atoms in optical lattices ul- 
timately led to a clean experimental realization of the 
superfluid-Mott transition [l|, ]2j]. This in turn inspired 
many intriguing developments marrying condensed mat- 
ter physics, quantum optics, and quantum information 
processing J3(]. Concurrent with the superfluid-Mott tran- 
sition, revivals of first-order phase coherence upon enter- 
ing the Mott phase in a 3D optical lattice have been ob- 
served |J|. In the ideal case of a perfectly homogeneous 
lattice, the revivals would be exact, i.e., reach the same 
coherence amplitude as the initial superfluid state, pro- 
vided the optical lattice wells were completely isolated 
from each other and from the environment (heat bath). 
The experiment [4j] however measures a significant damp- 
ing of the revival amplitude (so that approximately four 
to five revivals can be observed clearly). Here, the revival 
time t rcv is defined to be the time at which the first re- 
vival maximum of the long-range coherence (a^t)^^)) 
occurs after the sweep to the Mott phase. There are 
basically four mechanisms which might induce the ob- 
served damping: i) a finite tunneling coupling J of the 
lattice wells, ii) a finite coupling of the wells to the envi- 
ronment, and finally inhomogeneities in the lattice, such 
as Hi) an external trap potential (inducing finite-size ef- 
fects) , and iv) variations of the lattice depth and hence 
on-site repulsion U from site to site (i.e., diagonal disor- 
der). Here, we investigate the question of whether and 
how a finite tunneling coupling J of lattice wells induces 
a damping of the revival signal, i.e., mechanism i). In 
contrast to the other effects, this mechanism is an in- 
trinsic property of the system and survives in the ho- 
mogeneous continuum limit. It is also theoretically well 
under control and of fundamental interest, e.g., regarding 
the basic problem of equilibration, that is the question 
of whether and how fast the system approaches a local 
equilibrium state after the quench to the Mott phase. We 
point out, as these arguments already indicate, that the 
nature of revivals and their damping in an optical lattice 
is fundamentally different from the single- well case where 
the particle-number dependence of the chemical potential 
plays the major role [5j. Our strategy is to calculate per- 
turbatively the solution of the many-body problem on a 



general lattice, when the parameter 2irJ/U is small but 
finite, starting from the exact number basis solution at 
J = 0. Here J gives a typical scale of the hopping am- 
plitude, and U is proportional to the contact interaction 
coupling. We note that within the presently employed 
pertubative approach, we cannot discuss the late-time 
behavior but can extract the decay of first-order coher- 
ence from one given to the next revival cycle. 

Tunneling-induced damping of revivals. The Hamilto- 
nian we employ is of the conventional single-band Bose- 
Hubbard type with on-site interactions 



H = JM^al^ + -(al) 2 al, 



(1) 



where of, and a v are bosonic creation and annihilation 
operators at lattice sites fj, and v, respectively. We use 
the convention that summation is implicit over the equal 
Greek lattice site indices designated with an overbar. 
The matrix M^ describes the (possibly anisotropic) tun- 
neling rates on an arbitrary lattice [6j], and U is the on- 
site repulsion. The full operator equation of motion in 
the final state with given J then reads (ft = 1) 



. dcifj 
'~dt 



= JMnyCiy + UfluCln 



(2) 



where n^ — aj^a^ counts the number of particles per site. 
The exact solution of the above operator equation, for 
decoupled wells, J — 0, is given by 



a£(i) = exp{-it/n^K(0) 



(3) 



The revival time, where the argument of the exponential 
operator assumes integer multiples of 2ir, is t rcv = 2tt/U 
(for J = 0) . These coherence revivals exist because of the 
discrete spectrum of the filling operator ft M , i.e., because 
of the existence of particles, the "granularity" of matter. 
Calculations for collapse and revival of first- and higher- 
order correlation functions for the decoupled-well case 
J = were performed in [8{]. We mention here that in 
the opposite limit of hard-core bosons (U — > oo), hopping 
in an optical lattice from site to site, phase coherence 
collapse and revival may occur when an optical super- 
lattice is switched on [fj]. 



We now proceed to calculate the perturbative solution 
of |(2]) for J/U small but finite. Defining 



a M (i) = exp{-ii7n /i t}A /i (i), 



(4) 



we extract the dominant (for J/U <C 1) time-dependence 
from the full field operator a M (t), so that Ap(0) = 6^(0) 
and Ap(t rcv ) = a^(trev). The equation of motion 



id t Ap = JM mP exp{iC/(n M - np)t}A p 



(5) 



can be solved via an expansion into powers of the small 
parameter J (analogous to response theory with a per- 
turbation H int = JMppan&p), where the second-order 
solution reads 

A^t) = i M (o) - i [ dt'jM^^-^'A^o) 




t ft' 

dt' / dt"j 2 M^Mp P e lU{il ^-^ )t ' e 111 ^-^ 1 " A- P {Q) 
Jo 

+ 0(J 3 ). (6) 

Inserting the above second-order solution into the first- 
order correlation function evaluated at time t rcv yields 



.i;,(/„, ).u/,, v r) = (.i;,ii)LU0))-(^] x 



x {<Jit (0)<S(n M - n p )5(n ff - fip)Ap{0)^ M^M- pa 
+ 2 (^(°)<^ - fip)5(n a - hp)A p (0)} M 9p M pa 



+ 2 \ A p(0)S(n M - h P )5(m - h P )A a (0)j M^Mvp 

+ 0(J 3 ), (7) 

where 8{h v — n p ) is to be understood as an operator 
Kronecker delta, and the average is taken with respect 
to the initial many-body quantum state. The above ex- 
pression represents the full quantum result to second or- 
der in J, which is general insofar, within the single-band 
Bose-Hubbard model. This should be contrasted with 
the approach of [llj], which considers quantum collapse 
and revival in a two-dimensional (2D) respectively three- 
dimensional (3D) optical lattice, using resonant tunnel- 
ing in the reduced dynamics of a two-mode (2D) or three- 
mode (3D) model. 

The terms linear in J generally vanish in the final re- 
sult for the correlator ([7]), which is to be expected be- 
cause (A'(0)A a (0)) assumes its maximum value (= n in 
the homogeneous case) in a coherent state. Within linear 
response, the tunneling-induced damping therefore gen- 
erally vanishes. We now use such a coherent state, more 
precisely a product of coherent states at each site, which 
has Poissonian number statistics, to evaluate (7|), 



|coh) = niaV = Il' 



KI/2 




(8) 



where |n M ) are local number eigenstates, 
\np)rip and dp\a) 



\a) p a p with Poissonian distribu- 



tion and average (n^) = |a M | 2 . The state iJH) is the factor- 
ized many-body state after a quasi-instantaneous sweep 
starting deep in the superfluid state [6{]. We assume ho- 



mogeneity (no external trapping), a M = a, 



n, and 



calculate the two different terms occurring in the corre- 
lator (7]). For the relevant case of long-range coherence 
(i.e., the two lattice sites fi and a are far apart and do 
not share neighbors, the result then being independent 
of the distance of the two sites), the first term gives 



Al(0)S(hp - n 9 )5(n a - h p )A p (0)) M^M pa 



= « XXJ D 2 = ne-^ll{2n)D\ (9) 
Vfc=o / 

where D is the number of neighboring sites defined by 
D = J2v M^v for a fixed site /u, with the i/-sum running 
in the simplest case over only nearest neighbors to site /x 
(see however, discussion below of including next-nearest 
neighbors as well), and Iq is modified Bessel function; the 
Poisson probabilities to find k particles in a coherent state 
with the number average n are denoted Pk, n — e~ n n k /k\. 
Similarly (for long-range coherence), the last two terms 
in (7j) each give a contribution 

- (Al(0)5(n p - hp)S(n a ~ n p )A p (0)) M pp M- pa 



(OO \ OO 

5> 2 fc J D{D-l) + l -nY,Pl, n D 
fc=0 / fc=0 

= ^e-^I 2 (2n)D(D 1) + ^ e - Sn F {l<1} (n 3 )D , (10) 

where -Ffi.i}^ 3 ) is a hypergeometric function. The sec- 
ond term is generally negligible for D ^> 1, cf. Fig.Q] be- 
cause T,kLoPk,n ~ (EfctoPiL) 2 for the experimentally 
relevant range of n = 1 ... 10, with the ratio of J2"kLo Pk n 
t° EfcloPfe n) 2 slowly decreasing with increasing n. 

Collecting the correlation functions from © and (fit))) , 
the total result for the decay of first order coherence at 
the instant of the first revival reads (see also Fig. [I]) 



A(t rev ) 



4(0)i CT (0))-(it(t rev )i CT (frev) 



n 



= (^pj [D(2D-l)e-^I 2 (2n) +D e -^F {hl} (n 3 )] . 

(11) 

This represents our major result for a coherent state with 
Poissonian statistics, i.e., after a sudden sweep from deep 
within the superfluid phase, for which the original ground 
state ((HJ has no time to adjust. 

Numerical estimates and discussion. We now pro- 
vide a numerical estimate of l|lip for the experiment 
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FIG. 1: The result ifTTjl for the damping rate as a function of 
the filling n (bold line), for (2ttJ/U) 2 = 1.15X 10~ 3 and D = 6 
[j|. Dotted and dashed lines represent the contributions of the 
Iq and F{i,i} terms in ifTTjl. respectively. 



Q, which was performed in a cubic 3D optical lat- 
tice, with an initial depth Va = 8E r (corresponding 
to the superfluid phase), and a subsequent "jump" to 
a final depth Vb — 22E r (deep into Mott phase). In 
terms of the recoil energy E r = 2ir 2 / (mX 2 ) determined 
by the laser wavelength A, the tunneling coupling J in 
sufficiently deep optical latt ices c an be estimated by 
J/E r = 8{V/E r ) 3 / 4 exp[-2y/V/ET]/y^. Similarly, the 
on-site repulsion scales as U/E r = 4:\/2ir(V/E r ) 3 > 4 a s /\, 
where a s is the s-wave scattering length [7j . One obtains 
for the perturbative parameter of our calculation 



2vrJ 

IF 



8A 2 

^-exp 
ai 



(12) 



[4 



Assuming, like for 87 Rb, a s — 6nm, and a laser wave- 
length of A = 852 nm, we have for V = V B = 22E r the 
value {2nJ/U) 2 = 1.15 x 10~ 3 . With the number D = 6 
of nearest neighbors in a 3D cubic lattice and the experi- 
mental value [4j of the average filling, n = 2.5, we obtain 
A(i rev ) = 2.8 x 10~ 3 . Thus damping induced by tun- 
neling to nearest neighbors, evaluated on the basis of a 
superfluid state with Poissonian statistics cannot directly 
explain the amount of damping observed in experiment 
which is one to two orders of magnitude larger [lfj| |. 
The result (fTT)) however represents a lower bound 
for tunneling-induced damping. We now discuss vari- 
ous scenarios which would increase A(i rcv ) above this 
lower bound. Firstly, (fTTj) is the result obtained us- 
ing the Poissonian superfluid state (J> U) as the ini- 
tial unperturbed (zeroth-order) state. Assuming sub- 
Poissonian number statistics with a reduced number vari- 
ance increases the expectation values of Kronecker sym- 
bols S(h^ — h v ). In view of (V>| S(h^ — n v ) \ip) < 1 for 
all states \ip), we may obtain a crude upper estimate by 
just omitting all Kronecker symbols 8{n t _ l — n v ) in the 
expectation values. Inserting this upper bound, reason- 
able agreement with experiment can be reached, noting 
that EfcoPM.5 - 0-034 and (E^o^.s) 2 - 0.0385, 
which would then essentially be replaced by unity. Note, 
however, that actually realizing this upper bound would 



imply that the many-particle state approaches a product 
of Fock states with exactly n particles at each site - for 
which the coherence signal (A'„A a ) vanishes. Neverthe- 
less, a state "in-between" the coherent state and the Mott 
state could still display non-vanishing long-range coher- 
ence (though below its maximum), n > (A* A a ) > 0, 
whose revivals are more strongly damped than in the fully 
coherent case ifTTj). Actually, such an intermediate state 
with sub-Poissonian number statistics is automatically 
created by a superfluid-Mott quench with a finite sweep 
rate J6(]. The character of the final state (i.e., whether it 
is closer to the initial superfluid phase or the Mott state) 
depends on the time scale of the quench in comparison 
with the (inverse) chemical potential, i.e., the relevant in- 
ternal energy scale [a]. From the revival time of 550 usee 
[J|, one obtains for n = 2.5 an inverse chemical potential 
of /i _1 = (Un)^ 1 = 40 fisec. According to [4J|, the system 
is quenched in 50^isec from superfluid to Mott. The ini- 
tial V /E r = 8 and final values V /E r = 22 imply that 
J is quenched across approximately 3.7 e-foldings. In or- 
der to compare the results with the calculation in [0], we 
assume an exponential decay J(t) ex exp{— jt} for sim- 
plicity. This implies that I/7 ~ 13.5 /Ltsec, which is about 
a factor of three smaller than the inverse chemical poten- 
tial of [i^ 1 = (Un)^ 1 = 40 [isec. Therefore, the sweep 
is rather non-adiabatic, which corresponds to a small 
adiabaticity parameter v = Un/j sw 1/3 introduced in 
[Qj. However, the sweep is not sudden - for a sudden 
quench, the number fluctuations would essentially retain 
their initial value, (n 2 ) — (n) 2 = (Sn 2 ) = n for a coher- 
ent state. Due to the finite sweep rate J(t) ex exp{— jt} 
with v — Un/j ~ 1/3, though, the number variance is 
reduced by 60%, i.e., (Sn 2 ) w OAn. Even though this re- 
sult, derived in [6J applies, strictly speaking, in the limit 
of large fillings n ^> 1 only, one would expect a similar 
sub-Poissonian statistics also for n = 2.5, which would 
then increase the decay A(i rcv ) in ifTTj) significantly. An 
additional source for sub-Poissonian statistics could be 
the initial state itself, if we do not start deep in the su- 
perfluid phase. Closer to the transition line (but still 
on the superfluid side), the number variance decreases, 
(Sn 2 ) < n already in the initial ground state. (In the 
experiment [4], merely 60% of the atoms occupied the 
coherent condensate state initially.) 

Another effect which potentially increases A(i rov ) is 
to include couplings to next nearest neighbors, i.e., to 
effectively increase D. For the parameters of the ex- 
periment Q, the final next-nearest-neighbor couplings 
may be estimated as follows: Within a harmonic ap- 
proximation, the tunneling matrix elements are given 
by JM a/3 ~ J\ a -f3\ exp{-Tr 2 y/V /E r (a ~ /3) 2 /4}, where 
\a — P\ is the distance of lattice points a and (3 in units 
of the lattice spacing and J\ a -p\ depends polynomially 
on J a — p\ [12J, ll3j. There are 20 next-nearest (diagonal) 
neighbors in a 3D cubic lattice in addition to the 6 nearest 
(straight) neighbors. Due to the exponential reduction of 
their individual contributions when quenching deep into 
Mott phase (Vo/E r = 22), their contribution gives about 



a factor of two for the damping A(i rov ). On the other 
hand, using only a moderately smaller final Vq/ E r may 
make their total contribution to revival damping larger 
than that of the nearest neighbors. 

Concluding remarks. In conclusion, we have derived 
a rigorous second-order result ([7]) for tunneling-induced 
damping of phase coherence revivals in optical lattices 
sufficiently deep inside the localized Mott phase 2ir J <C 
U. Evaluating expression (7|) for a coherent state, we 
obtained a lower bound for tunneling-induced damping 
(fll"]) . which is too small to explain the experiment Q. 
However, the incorporation of next-nearest (diagonal) 
neighbors as well as sub-Poissonian number statistics in- 
duced by the initial state and the finite sweep rate signif- 
icantly enhance the tunneling-induced damping A(t rBV ) 
such that - even though it does perhaps not fully repro- 
duce the experiment [4j - it should constitute an observ- 
able fraction of the measured decay. Unfortunately, the 
precise value of the tunneling-induced damping of phase 
coherence cannot be derived from the information given 
in [J] since part of the relevant input - such as the exact 
dynamics J(t) - is missing. Thus, while the dephasing 
mechanisms i) to iv) discussed in the Introduction are 
qualitatively well understood, it is not possible to pre- 
cisely disentangle their quantitative contributions from 
the data at hand. 

While our prediction cannot be compared accurately to 
existing experimental results, a relatively modest modifi- 
cation of the parameters in an experiment like that of Q 
will allow for its test. In particular, the exponential de- 



pendence of the tunneling-induced damping A(i rov ) oc J 2 
on the laser amplitude y / Vo/-E'r 1 _implying 7.4 e-foldings 
for (2ttJ/U) 2 in the experiment [4|, distinguishes (at low 
temperatures) this mechanism i) from other sources like 
inhomogeneities due to external trapping Hi). 

We finally point out the importance of tunneling cou- 
pling for the equilibration of many-body states on the lat- 
tice. Nonequilibrium states of the Bose-Hubbard model 
were studied for various cases, e.g ., quenching from the 
superfluid to the Mott side [fj, ]l5j]; from Mott to super- 
fluid [la ]: or for hard-core bosons in superlattices [9|]. 
While our approach is perturbative, yielding first-order 
coherence up to second order in 2ttJ/U, it provides a 
first estimate on how important coupling of the wells can 
turn out to be for an equilibrium or non-equilibrium state 
to be established after a quench. Taking the decay of 
the first-order coherence as an indicator for locally ap- 
proaching the equilibrium state, we conclude that equili- 
bration occurs faster for many neighbors and for a larger 
tunneling rate (as one would expect) but also for sub- 
Poissonian number variance - this is somewhat surpris- 
ing, as one would generally expect that a state which is 
already closer to the Mott state decays slower. 
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